Mass distribution from quark matter equation of state 
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We analyze the equation of state in terms of quasiparticles with continuously distributed mass. 
We seek for a description of the entire pressure - temperature curve at vanishing chemical potential 
in terms of a temperature independent mass distribution. We point out properties indicating a mass 
gap in this distribution, conjectured to be related to confinement. 
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According to the proposal of A. JafFe and E. Wit- 
ten, a successful quantum Yang-Mills theory must have 
a mass gap In heavy ion collisions a deconfined 

phase is expected to form, and the produced quark-gluon 
plasma (QGP) is described by Quantum Chromodynam- 
ics (QCD). Thus one can expect the appearance of such 
a mass-gap in the spectral function of the basic QGP 
degrees of freedom, namely quarks and giuons. In this 
paper we perform a quantitative analysis on some results 
from lattice QCD on the equation of state (eos) and re- 
construct it from a mass distribution of non-interacting 
quasi-particles. We present strong indications for a mass 
gap in this distribution. 

We have used earlier a mass distribution for massive 
quarks and developed a coalescence picture 3 to de- 
scribe hadronization of deconfined quark matter, and re- 
produced final state hadron ratios and transverse spectra 
successfully. That model was based on an earlier coales- 
cence model 0, 3], where quarks and giuons had finite 
effective masses without any width. The apparent en- 
tropy reduction problem by coalescence with an associ- 
ated reduction (confinement) of color degrees of freedom 
can be resolved by assuming sufficiently massive partons 
around the hadronization temperature in the precursor 
matter. The necessary mass scale for quarks is about 
300-350 MeV and even higher (about 700 MeV) for giu- 
ons 0, thus we could assume that in the prehadroniza- 
tion stage the heavy giuons decay into quark - antiquark 
pairs [3|. Recently, partonic level models of heavy ion 
reactions also utilized the quark coalescence picture suc- 
cessfully 0,0]. 

Considering quark coalescence as the mechanism of 
hadronization one has to deal with the question, how 
to make a hadron with a mass lower than the sum of two 
parton masses. In order to solve this problem we have in- 
troduced distributed mass partons into our hadronization 
model 01 ■ Having in medium partons in quark matter as 
precursors of emerging hadrons in mind, we connect now 
the distributed mass parton picture to a simplified treat- 
ment of spectral functions. 



The equation of state of an interacting system, when 
analyzed in terms of quasiparticles, is coded in the spec- 
tral function p{uj,p). Our ansatz to this assumes a par- 
ticular form: 

p(c.,pl = 27r^e(m2)(e(c.)-e(-c.)) (1) 

with m? = up' — 'j? and Q(x) being the step function. 
A continuous w{m) mass distribution describes a finite 
width ansatz for the spectral function. The normaliza- 
tion of the spectral function, p{uj,'p)Loduj /it = 1, re- 
quires w{m)dm = 1. The conventional quasiparti- 
cle approach on the other hand often explores a Breit- 
Wigner form of the spectral density [1], 
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{lu - EY +72 {uj + EY + 72 



(2) 



with E^ — A/2 — 72 -|- and temperature dependent 
parameters M(T) and 7(T'). 

Thermodynamical consistency of the quasiparticle pic- 
ture imposes constraints on the mass distribution, w(m), 
in particular on its dependence on the temperature or on 
other medium parameters In this letter we investi- 
gate the possibility of a temperature independent mass 
distribution and therefore neglect the mean field term 
for consistency. The total pressure at vanishing chemical 
potential is given as the following integral: 



poo 

p{T) = / w{m) p{m,T)dm. 
Jo 



(3) 



One may suppose that only a single mass scale occurs in 
the mass distribution, so it can be expressed by a dimen- 
sionless distribution: 



(4) 



The normalization integral for w is inherited by the shape 
(form factor) function f{t): 



w{m)dm = / f{t)dt = 1. 
Jo 



(5) 



The quark gluon plasma at vanishing chemical poten- 
tial has the pressure 
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p{T) = a{z)KT\ 



(6) 
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with z = Tc/T and k being the Stefan-Boltzmann 
constant. In the Boltzmann approximation the fixed 
m-contributions are given by the Bessel K- function, 
p(m,T) oc T^<^>{m/T) withll^ ^{u) = u'^K2{u)/2. De- 
viations in ^{m/T) ~ p{m, T)/p{0, T) due to using Bose 
or Fermi distributions as a function of m/T never exceed 
six per cent. Thus in the distributed mass model the <t{z) 
function in this approximation is given by the integral 



a{z) 



' fit)^-^K,izt)dt. (7) 



This may be recognized as the so called Meijcr 
K-transform [ll|, [l^ (a generalized Laplace transform) 
of the f{t) function. The inverse of this transformation 
yields the mass distribution function in terms of the ob- 
served a{z) values : 
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'^'°°^(,)Mf!)dz. (8) 



This raises a peculiar question: is it possible to show, 
that to any a{z) function extracted from an equation of 
state (e.g. from lattice QCD calculations) there exists 
a unique mass distribution f{t) with the mass scale pa- 
rameter kept temperature and chemical potential inde- 
pendent? In this case the shape of the mass distribution 
is not arbitrary. Of course, the Meijer K-transform is in- 
vcrtiblc, but one has to check whether the f{t) function 
obtained by eq. jS]) is positive semidcfinit and normalized 
to unity. The normalization is the easier problem, the 
ct(0) limit being directly the integral of the f{t) func- 
tion due to the small argument behavior of the Bessel 
K-function. It can, however, be difficult to arrive at 
a nowhere negative f(t) by knowing cr(z) only at some 
points on the real z-axis. 

Before investigating any particular ansatz for a{z) let 
us consider an important general property. There is a re- 
lation between the integration moments of this quantity 
(the scaled pressure) and the normalized mass distribu- 
tion, f{t): 

pQO poo 

M„ = / z"-V(z)dz = /„ / f{t)t-"dt (9) 
Jo Jo 

with 

In = \j^ u-+'K,iu)du = 2'T (2 + ^) r (10) 

where T{x) is Eulcr's Gamma function. This is finite for 
positive n and divergent for zero or negative integer val- 
ues. We conclude that as long as the M„ moments of 
the eos curve are finite so must be the inverse mass mo- 
ments of the mass distribution. Since due to construction 
ct(0) = 1 and a{z) is rapidly decreasing due to confine- 
ment for large z = Tc/T (low temperature), any mass 
distribution reconstructing the equation of state of QCD 
must be suppressed for low masses. 
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FIG. 1: The pressure normalized to the massless Stefan- 
Boltzmann value as a function of the scaled temperature T /Tc 
for different constant mass relativistic gases (full lines in order 
for M/Tc = 0, 1, 2, 3, 4, 5, 6 and 7) and from lattice QCD data 
of Ref.[l3| (down triangles) , of Ref. [TJ] (open circles) and of 
Ref. [iq (up triangles). 



Let us now discuss how to obtain a particular func- 
tional form for a{z). The high-temperature (small z) 
expansion of $(zt) leads to 
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(11) 



with 7 being the Euler-Maschcroni constant. Acciden- 
tally the perturbative QCD pressure (taken as the finite 
part at the scale 27rT) shows a similar structure, 



V{T) 



1 - a2g^ -I- a^g^ + h^g^ In 



27rr 



(12) 



kT^ - . . - ^ 

(02 « 0.072,04 « 0.061,64 sa 0.008 for A^/ = 3 based 
on Ref. [l3|). It is possible to fit this form as a high-T 
asymptotics by assuming a scaling of expectation values, 
like (m2) = cg^T^, (m^) = c'^'^T^, etc. This is the basis 
of the traditional quasiparticle picture [lB|, at the same 
time it predicts a width changing with the temperature. 
This view assumes a temperature-dependent mass distri- 
bution, w{m, T), which - for the sake of thermodynami- 
cal consistency - would require a temperature-dependent 
mean field pressure, —B(T) to be taken into account. 

There is, however, another possibility, which we would 
like to pursue in the present article. The low-argument 
expansion eq. (jlip fails if the expectation values, like 
(jTi^), (m**), etc. are divergent. In fact this assumes 
a high-mass tail of the wirn) distribution not decaying 
faster than m^^ . As we shall point out later, our nu- 
merical efforts to obtain w{m) agree with this statement. 

Fig[T] presents the normalized pressure for relativistic 
Boltzmann gases with several fixed masses. The lattice 
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QCD eos data of the Budapest- Wuppertal group [l3| and 
of the Bielefeld group IJ] seem to lie everywhere below 
the curve for the mass M = 3Tc, recent MILC data (isj 
below the curve for M = 2.5Tc. As a consequence, if one 
accepts this property also for lower temperatures where 
actually no reliable simulations are available, the mass 
spectrum w{m) would not contain any mass lower than 
M = 3Tc or M ~ 2.5Tc, respectively. This property can 
be verified by rigorous mathematical estimates of upper 
bounds for w{m) on the interval < m < M j2^ . 

The pressure of hot QCD has been recently calculated 
up to 0{g^ln{l/g)) [17|. The result contains formally 
ln(27rT/A) terms, but according to the suggestion of the 
authors the coupHng g{A) should be taken at A « 6.47T 
[l^ and this way the temperature dependence of the nor- 
malized pressure stems from the temperature dependence 
of the coupling strength, renormalized relative to a scale 
proportional to the temperature. Agreement with lat- 
tice QCD eos data is achieved at the highest computable 
level only, with an extra fit of a constant which is not 
calculable perturbatively. 

It is intriguing that for practical purposes the 0{g^) 
formula by using the 1-loop renormalized g{T) ~ 
l/61n(T/A) can be also fitted to lattice data by fitting 
A. This results in a formula 



PjT) 
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K 



(13) 



reaching zero pressure at T « T^. Fits to different lat- 
tice QCD equations of state leads to quantitative, but no 
qualitative differences. For the data of Ref. [l^l we ob- 
tain K = 0.54^7? = 1.76, for Ref.[l3 K = 0.43, r? = 1.6 
and for Ref.^ K = 0.22, r; = 1.12 (cf. dotted lines on 
FigEl). 

There are some theoretical signs on the other hand, 
that the g{T) = l/&ln(T/A) formula is not necessar- 
ily the really high temperature limit prediction of QCD. 
The finite temperature renormalization group result is of 
type0 



a(Q^r2) 



y2 



^ - ^ ) (14) 



with b = l/a{Ql,T^) + bo and bo ^ (llA^c/2 - 
2Nf /3)/{4:Tt) being the one-loop perturbative beta func- 
tion coefficient. This is usually considered in the Qq ^ 
limit and then, assuming a sharp thermal distribution 
of values, ~ [aixT)"^ is taken. This leads to the 
conjecture 
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a((a7rr)2,r2) ^(y2) 



bo In 
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(15) 



Since Qo was large and A is around Tc, also the coeffi- 
cient a is taken as a large number. None of these assump- 
tions is established by the QCD itself. The assumption 
Qo ^ contradicts the T oo limit, the spread of a 



thermal distribution of possible values also increases 
like T^, and finally there is always a non- negligible in- 
fluence of low physics on the coupling at any flnite 
temperature. In fact calculating the thermal distribution 
of /T^ between two massless, Boltzmann-distributed 
particles one obtains easily 
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for > 0. The probability of having = is finite at 
any temperature, P(0) = 3/64 and this is the maximum 
of P{x). As a consequence higher twist effects which are 
not infrared safe (like InQ^, etc.) might destroy the 
often quoted logarithmic scaling of the coupling constant 
with temperature. It is probably the best to consider 
a K / lTi{r]T /Tc)-type formula fca fT5)) as a standard, but 
not the only possible fit to the lattice eos, with some 
parameters of nonperturbative origin. 

In order to evaluate the inverse Meijer transform, 
eq.®, one has to approximate the lattice QCD data by 
an analytic 17(2;) function. A family of mass distributions 
can be Meijer-transformcd analytically: 
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a{T) 
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A^)' "e(m-A), 



r(^) 



(17) 



with < < 2. The ly = 2 limit belongs to a Dirac-delta 
mass distribution, the 1^ = limit to the Bessel function 
A'o with logarithmic asymptotics for high temperature 
(small argument). All these ansatze contain a mass-gap, 
the distributions being zero for m < X. The v = 1/2 
value leads to the particularly simple eos: a = e~^^^ . 

We found that an overall fit in the range of known 
lattice data is also achieved by the analytic ansatz 



(t(z) = cxp (— Az) 
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1 + ef^^")/*" 



(18) 



with z = Tc/T, A = 1.05, a = 0.90,6 = 0.11 for data 
from [H, A = 0.87, a ^ 0.90, 6 0.10 for data from [Til 
and A = 0.56, a = = 0.10 for data from We 

note that in Refs.[ll, [li] « 170 MeV, but in ReL^\ 
Tc ~ 190 MeV was taken. These fits are demonstrated 
in Figl2] where the different sets of lattice QCD data are 
compared with the fitted (j{z) = p/psB curves of eq. lfTS]) . 

In this letter we investigate the lattice QCD eos data 
of Rcfs. [Tsj l closely, but they qualitatively agree 

with other results on this issue. The rise at moderately 
high temperatures (low z) cannot be accommodated by 
quantum statistical effects, but it can be characterized 
as the effect of an exponential factor exp{—XTc/T) in 
the range from to 2.5Tc (cf. FigEJ. While this mod- 
erately high-temperature behavior is well fitted by the 
pure exponential (7{z) function, the part below Tc is more 
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FIG. 2: The lattice QCD pressure normalized to the massless 
Stefan-Boltzmann value as a function of the temperature T 
from lattice eos results of Ref. [H (above), Ref.fTj] (middle) 
and Ref. [3| (bottom). Our fits are indicated by the continu- 
ous lines, the 1 — K /ln{riT /Tc)-type fits by the dotted lines. 



reduced. The <j{z) = 1 — K/ln{rj/ z) form is also able to 
fit T > Tc data, but it goes to negative values at a finite 
temperature, which is unphysical. Our exponential fit is 
overall positive. A most satisfying extrapolation would 
interpolate between these two functions. 

In FigH lattice data from Ref.[ll| (a), Ref.[li| (b) and 
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FIG. 3: The mass distribution function obtained by eval- 
uating the complex integral from eq.© (boxes) and using 
the analytic fit ea. (|18|l to lattice QCD eos data of Ref. [l3| . 
The full line corresponds to ea. (|20|l . circles show the curve 
obtained by using the 1 — K /ln{iiT /Tc) type fit. 



Ref. (c) on p/psB = (^{z) f-s a function of the temper- 
ature T /Tc are plotted. In a one-loop resummed pQCD 
motivated approach, using a mass directly proportional 
to the temperature the approach to one is logarithmic, 
1 — K/ \n{riT/Tc) (dotted lines). The exponential behav- 
ior, on the other hand, supports the presence of a lowest 
mass in high-temperature QCD. While this fact depends 
on the low-tcmpcrature drop of the pressure curve, it is 
not easy to consolidate the effect due to quantitatively 
different pressure curves presented by different lattice 
QCD calculations. Although we do not intend to re- 
view lattice QCD eos calculations in this paper, we note 
that the investigated simulations differ in the correspond- 
ing value of the physical pion mass (m^ « 140 MeV for 
13 1, « 540 MeV for Q and « 300 MeV for 
15|). There can be further differences of technical na- 
ture, which we do not feel to be able to comment on. In 
our further analysis we choose the data of the Budapest- 
Wuppcrtal group [isj to seek for a corresponding mass 
distribution, but of course the same exercise can be done 
for other sets of pressure data, too. 

To evaluate the integral given by eq.dH), we choose a 
simple path parallel to the imaginary z axis, z = c + ito. 
With numerical integration we obtain an f{t) mass dis- 
tribution shown in FigO by full boxes. Fluctuations at 
small masses are due to limitations of the applied numeri- 
cal method. The part of the mass distribution shown here 
reconstructs the T > Tc part of the pressure curve nicely, 
but it fails to approximate the pressure aX T < Tc- In 
the following we seek to understand this phenomenon. 

One can obtain simple analytic approximations for 
the f{t) function by expanding the expression for cr(z), 
eg. psp . However, requiring a convergent expansion, one 
arrives at two distinct scries expansions: one for z < a 
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and another one for z > a. 

e^""^ + . . . for z < a 



a{z) 



for z > a. (19) 



The inverse Meijer K-transform of the simple exponen- 
tial, exp(— Az), can be given based on an analytically 
known integral (cf. eqlT7]for v = 1/2): 



(20) 



The above expression is valid for t > A, for smaller 
t ~ m/Tc values f{t) is identically zero. Hence the t- 
integration in the Meijer K-transform, when determining 
the pressure contribution, starts at t = A. Physically 
this corresponds to a lowest mass in the continuous spec- 
trum, to a mass gap. Since both the approximations to 
T < Tc and to T > parts of the pressure contain a 
leading exponential factor (A and A -I- 1/& respectively, 
cf. ea. (fT9|) ). eos data seem to support a lowest value of 
a continuous mass spectrum both in moderately low and 
moderately high temperature quark matter (see the full 
line in FigE]). 

Actually, requiring z > a is equivalent to a Hagedorn 
limiting temperature Th = Tc/a, and in fact transforms 
back nearly to an exponentially rising mass spectrum 
part. In this regime the QCD matter also has been fitted 
by a hadron resonance gas [l^ . 




FIG. 4: The normalized lattice QCD pressure and the pres- 
sure fitted to convergent series expansions, eq. (|19|l . obtained 
by numerically re-integrating f{t) functions given by ea. (|20|l . 

Substituting the respective /(t)-s for T < Tc and 
T > Tc into ea. (|20|) we calculate the pressure from eq.([7]). 
These two curves are shown in Fig|4l together with the 
lattice QCD results of Ref.[l^. A numerical method de- 
signed to obtain an overall non-negative (probability like) 
f(t) distribution, which fits well some a{zk) = Sk points, 
is represented by the maximum entropy method (MEM). 



We applied this method to the lattice QCD eos data dis- 
cussed in this paper in order to obtain a mass distribu- 
tion: both by using a MEM program designed to invert 
the Meijer K-transform and also by searching numeri- 
cally for the inverse Laplace transform of a{z) first. We 
failed, however, to obtain better numerical results then 
by evaluating the complex integral eq.® as discussed 
above. 

For calculating quark number susceptibilities or higher 
order Taylor coefhcients the Boltzmann approximation 
becomes unreliable; starting at the fourth order the 
Boltzmannian term no longer dominates the Fermi dis- 
tribution. Experience with the hadronic resonance gas 
model supports the expectation that dependence on 
chemical potential also can be interpreted in terms of 
clustered but non- interacting components 1^, 2^, 21 1 . 

In conclusion, we have analyzed lattice QCD pressure 
data in terms of a continuous, temperature independent 
mass distribution. We find a strong indication for a finite 
mass gap in such quasiparticle models, the details de- 
pending on the low temperature behavior of the pressure 
curve. Since all simulation data are below the AI ~ 2.5Tc 
curve, the immediate conclusion would be that the p{T) 
curve can be fitted by components with higher mass only. 
Allowing for a milder drop of the pressure at low tem- 
perature the lowest mass may be lower, we presented 
an example with Al « XTc with fitted A-values near to 
one. In general for any p{T) curve showing finite T"'+'^- 
weighted integrals for p/po the low-m behavior of w{m) is 
restricted by finite integrals of mr'^w{m). Since a single- 
mass p{T)/pq curve cannot fit the lattice QCD equation 
of state obtained by any of the groups calculating it, these 
data demand a finite width mass distribution. 

For the physical problem of quark matter we have 
learned from the above analyses that either the mass dis- 
tribution is temperature dependent and then the thermo- 
dynamical description is rather complex then, or there 
is a mass gap compatible to the equation of state un- 
less the pressure rises again at low temperatures (where 
we have presently no simulation data, but the idea of a 
non-interacting pion gas would correspond to a pressure 
higher than zero). 
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